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Abstract 

Coincidences between charged particles emitted in the /3-decay of ^^Li were 
observed using highly segmented detectors. The breakup channels involving 
three particles were studied in full kinematics allowing for the reconstruction of 
the excitation energy of the ^^Be states participating in the decay. In particular, 
the contribution of a previously unobserved state at 16.3 MeV in ^^Be has been 
identified selecting the Q;+''He^ a+^He+n channel. The angular correlations 
between the a particle and the center of mass of the ^He+n system favors spin 
and parity assignment of 3/2^ for this state as well as for the previously known 
state at 18 MeV. 
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The study of the /3-decay of ^^Li constitutes a challenge for experimental- 
ists. The high value (20.557(4) MeV 0,0]), combined with the low particle 
separation energies in ^^Be, opens a plethora of decay channels, /3-7 /3-n 
/3-2n /3-3n /3-a 0, /3-t and /3-d These decay channels involve a 
variety of emitted particles over a wide energy range, which only recently has 
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been open for detailed spectroscopic studies due to the availability of segmented 
detectors. The /3-n channel, corresponding to the lowest energy threshold, at 
0.504(6) MeV, has recently been studied in a series of /3-7-n coincidence ex- 
periments and /3-delayed 7 doppler broadening studies 11, [l2|, which 



mapped the states of ^^Be fed in the ^^Li /3-decay up to 10.6 MeV excitation 
energy. Information about the spin and parity of some of these states has also 
been obtained [loj . The remaining /3-delayed channels include the fragmenta- 
tion of ^^Be states into more complex channels, in particular the five body 2a3n 
channel makes the identification of the original state in ^^Be difficult. 

The first /3-delayed charged particle emission study done in coincidences and 
with time of flight particle identification established the presence of two different 
decay channels, three-body na^He and five-body 2a3n [l3|. Two states in ^^Be 
were proposed from energy considerations to contribute to these decays, one 
at 10.59(5) MeV identified in reaction studies [3| and a previously unobserved 



state proposed at ~18.5 MeV [13|. The second state was later confirmed by an 



experiment using a gas-Si telescope for particle identification [15l |. The energy 
spectra of the ^'^°Be, "^'^He and d,t decay products were studied and the rec- 
ommended energy for the state was 18.15(15) MeV. This state will be referred 
from here on as the 18 MeV state. 

Furthermore, there have been experimental hints of /3~ feeding to states 
in between these two states in ^^Be. In the gas-Si telescope experiment it 
was reported [3| thai a peak at 1.2 MeV in the ^Be-f ^°Be spectrum could be 
interpreted either as neutron emission from a ^^Be 14.5 MeV state to the ground 
state in ^°Be, or as neutron emission from the ^^Be 18 MeV state to the first 
excited state in ^°Be. A later experiment [lH that studied the neutron-recoil 
doppler broadening of 7-lines following ^^Li decay reported a 6% of /3-feeding 
to excited states in ^^Be in the region between the 10.6 and 18 MeV states. 

Additional information on the ^^Be states above the charged particle thresh- 
olds has been obtained by two recent reaction studies 1^, 17 1. From the spin- 



parity assignment given to these new states they are not expected to be fed in 
the /3-decay of ^^Li. 

On the theoretical side, shell model calculations [3] predict that a significant 
amount of the /3-decay strength lays in the region between 10 and 18 MeV in 
^^Be. Certainly, the shape of the Gamow- Teller strength distribution depends 
on the s/p mixing of the ^^Li ground state wave function. 

These results underline the importance of studying the ^^Li /3-delayed charged 
particle emission in great detail. In this Letter we present the study of the ^^Be 
excited states, fed in the ^^Li /3-decay, using the 
-'^-'^Be*^a-|-^He^n-)-Q;-|-^IIe channel, previously proposed in Q. We have es- 
tablished the presence of this channel in the breakup of the ^^Be 18 MeV state 



13l . Il5l |. Complete kinematical analysis is possible from the detection of the 



energy and direction of the emitted particles provided by the use of segmented 
detectors, which allows to reconstruct the ^He-|-a channel to obtain the exci- 
tation energy spectrum of ^^Be. Here we find the contribution of a previously 
unidentified state at 16 MeV, and determined, from the angular correlation 
study, the spin and parity of both states. 
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The experiment was carried out at the ISOLDE PSB FaciUty at CERN, 
where the ^^Li activity was produced in fragmentation reactions in a Ta target 
irradiated with 1.4 GeV protons from the CERN PS Booster. The target con- 
tainer was connected to a surface ionization source and the produced ions were 
extracted with a voltage of 30 kV. The 1+ ^^Li ions were mass separated from 
the different isotopes using the general purpose separator (GPS). The ^^Li beam 
was subsequently brought to the center of the setup where it was stopped in a 
50 jig/cve? carbon foil. The use of a thin foil, which minimizes the energy loss 
of the emitted charged particles, was possible due to the low extraction volt- 
age used in the experiment. The production in our experimental chamber was 
estimated to be 800 ^^Li ions/s, obtained from the daughter ^^Be f3a activity 

The experimental setup consisted of three particle telescopes, with 60 /im 
thick Double Sided Silicon Strip Detectors (DSSSD's) backed by Si pads for (5 
detection. Each DSSSD has its charge collection divided in 16 vertical and 16 
horizontal strips, obtaining an angular resolution of ~3.5° per 3x3 mm^ pixel. 
The telescopes were mounted in fixed positions on the surface of an aluminum 
table. The two detectors facing each other covered solid angles of 4% and 6% 
of 47r, respectively [20|. 

The silicon detectors were calibrated using ^^^Gd and the standard triple 
alpha source (^'^^Pu, ^^^Am and ^^*Cm). An analysis method developed within 
our collaboration [2l| takes advantage of the highly segmented nature of the 
DSSSD's to precisely define the path of the detected particle and thus allowing 
for the reconstruction of the energy losses in the non-active layers. The combined 
effect of the dead layers and electronic noise limits the detection threshold to 
160 keV for alpha particles. 

We concentrate here on the charged particle coincidences recorded in op- 
posite detectors, setting the angle between the detected particles between 120° 
and 180°. This geometry is optimum for the study of three body processes 
by detection of at least two of the products in coincidences. In our case 90% 
of the coincidences were recorded in the pair of opposite detectors. The time 
window for the analysis was set to the first 60 ms after the proton impact, time 
long enough to allow for the diffusion of the products from the target. This 
time window enhances the observation of the ^"'^Li /3-decay over the decay of the 
descendants. The directions of the detected charged particles were determined 
from pixel identification, which allowed for the calculation of their momenta as 
well as the deposited energy. In case of a three body process one can recon- 
struct the energy of the third undetected particle, the neutron. Assuming that 
the particle of lowest energy is the most massive one, the individual neutron, 
a and ^He energies are plotted versus excitation energy in ^^Be in Fig. [H^a). 
The advantage of this scatter plot, for sequential breakup of broad states, is 
that the first emitted particle will be along a line whose slope and offset are 
given by the mass ratio of the fragment to ^^Be and the Q- value of the breakup, 
respectively 22| . In the plo t most of the events are a-a coincidences from the 
five body 2a3n channel 3) H]j which are not correctly assigned in this plot as 
the three neutrons are not detected. Therefore, the events associated with the 
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2a3n channel are distributed over the plot and act as background for the other 
channels. There are two distinct features in Fig. [TJa). First, a horizontal line 
at 10.6 MeV excitation energy, associated to the breakup of the 5/2^ 10.6 MeV 
state in ^^Be [l^. Second, a grouping along the hue of 11/7 slope and with 8.35 
MeV offset. In this type of scatter plot, in case of sequential breakup through 
an intermediate resonance of energy Qji and mass M/j, the energy of the first 
emitted particle Ei will appear along a diagonal line following E=Qij+ -jg^^Ei, 
where M and E are the mass and excitation energy of the decaying state. In 
our case Or corresponds to the ''He+a separation energy in ^^Be, Q_r=8.35(2) 
MeV 0, [23|. Therefore, the line in Fig. [Ija) corresponds to a particles emitted 
in the breakup of ^^Be states through the ground state of the ''He resonance 
There are two regions with accumulations of events on this line, around 16 
MeV and 18 MeV excitation energy, indicating the presence of states in ^^Be at 
these energies. 

We selected events in the ^^Be*— ^^IIe(gs)+a^^Hc+a+n channel by gating 
on the charged particle with the highest energy, which is assumed to be the a 
particle. The two dotted lines in Fig. [TJa) show the interval analyzed, between 
the lines of slope of 11/7 and offsets of 8.93 and 7.93 MeV. The "Be excitation 
energy spectrum for these events, obtained from the sum of n, a and ^He energy 
plus Qna^Hej is shown in Fig. [Tfb). The two groupings of events observed in 
the scatter plot of Fig. [IJa) correspond to the clear peak around 16 MeV, plus 
the contribution at higher energies of the previously observed -'^^Be state at 18 
MeV ^Ml 

The peak at around 16 MeV indicates the presence of a state in 



^^Be previously unreported in the literature. 

We performed a Monte-Carlo simulation of the breakup of the 16 and 18 
MeV states in ^^Be into the three body na^He and five body 2a3n channels, 
the latter channel considered as a broad background. The resonances were de- 
scribed using the single-channel non-interfering R-matrix formalism as described 
in the Appendix A of Ref . . The a particle penetrability was calculated us- 
ing the orbital angular momentum values compatible with the spin-parity of the 
"'^"'^Be states and the 3/2~ ground state in ''He. The penetrability of the neutron 
connecting the ^He ground state with the O"*" ground state in ^He was calcu- 
lated for 1=1 orbital angular momentum. The five body 2a3n breakup for the 
two states was modeled following the prescription given in Ref. Q. The sim- 
ulation including the contribution of the different channels is shown as the red 
histogram of Fig. W[h), with each component color-coded following the legend. 
The parameters used to describe the different resonances are listed in Table [1] 

The resonance parameters were evaluated via tests of the ^^Be excitation 
energy spectrum. We performed a series of Monte-Carlo simulations varying 
the energy centroid and the reduced width individually, fitting the resulting 
values to a parabola to obtain the minimum. We obtained Eo=16.3(l) MeV 
and 72=0.05(1) MeV, with a full width at half maximum (FWHM) of 0.7(1) 
MeV, for the first state and Eo=18.4(3) MeV and 7^=0.11(4), FWHM=1.6(6) 
MeV, for the second. For both resonances the channel radius was set to 4.9 
fm, r=ro(Ay'^+A2^'^) with ro=1.4 fm, and the resonance shift parameter was 
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calculated using the formulation of the single-channel non-interfering R-matrix 



formalism as described in the Appendix A of Ref. [2j] . The energy and width of 
the latter state were previously determined in [l5| from the fit of singles spectra. 
Values of Eo=18.2(2) MeV and FWHM=1. 25(15) MeV were obtained from the 
9'i°Be spectrum; Eo=18.0(l) MeV and FWHM=0.8(1) MeV from the ''■'^He 
spectrum and Eo=18. 15(15); and FWHM=1.3(3) from the triton spectrum. 
The resonance parameters obtained in this work coincide within error bars with 
the previously published parameters, although are more in line with the values 
obtained from the ^'^°Be spectrum. Analysis of the triton spectrum from the 
same experiment yields an energy and width of this state of 18.35(30) and 1.5(4) 
MeV respectively [l^], consistent with the results obtained from the na^He 
channel presented here. 

The fact that the three body breakup channel is sequential and occurs 
through an intermediate state, ''He, with J^0,l/2 permit the use of the angular 
correlations to deduce the spin of the initial state in ^^Be. The angle 9 between 
these two directions is distributed according to the expression l-|-A2(2cos^ 0- 
l)/2, where A2 depends on the spin of the initial and intermediate resonances. 
As the ^^Li ground state is 3/2~ [2^, the selection rules for allowed transitions 
indicate that mainly l/2~, 3/2~ and 5/2~ states in ^'^Be are populated, whereas 
the spin and parity of the intermediate ''He(gs) resonance is 3/2~ j2^. We cal- 
culated the appropriate A2 parameters from the parameterization given in [S^l , 
obtaining A2 = l, and -0.714 for 1/2, 3/2 and 5/2 respectively, and performed 
a test of the simulation of the angular distribution of the breakup assuming 
every possible spin. In order to separate the contribution of each individual 
level we gated the sum energy spectrum in the following way: (i) the 15.7 to 
16.6 MeV interval is chosen for the characterization of the 16 MeV state, as its 
contribution in this region is higher than the contribution of the 18 MeV state, 
see Fig. [TJb), (ii) the gate for the 18 MeV state was chosen between 17.4 and 
18.6 MeV, asymmetric with respect to the known resonance centroid due to the 
weight of the Fermi factor. The experimental angular distributions for both 
cases are shown in Figs. HJb) and (a) respectively, with the simulations for the 
most favored spin assignment overlaid, and color coded according to the legend. 
One has to keep in mind that the spin of the 18 MeV state had to be studied 
first, as part of the statistics in the region of the 16 MeV state comes from the 
tail of the 18 MeV state. The results of the test, for 19 degrees of freedom, 
for the 18 MeV state were 67.54 for 1/2", 16.62 for 3/2" and 30.62 for 
5/2~, which show that a spin of 3/2~ is favored. In the case of the 16 MeV 
state, the results of the tests, setting the spin of the 18 MeV state to 3/2~, 
again for 19 degrees of freedom, are 70.50 for 1/2^, 20.57 for 3/2^ and 
58.68 for 5/2~, leaving 3/2~ as the only possible spin. In summary, the analysis 
of angular correlations in the ^He(gs)-|-a channel supports an assignment of spin 
of 3/2^ for both states in ^^Be. 

The /3-decay branching ratios to the ^^Be states discussed in this Letter 
were obtained from the Monte-Carlo simulations, rather than directly from the 
coincidence data. This is due to the difficulty in determining the experimental 
coincidence efficiency for the different channels discussed, as it depends on the 
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kinematics of the channel in question. The Monte-Carlo simulation, including 
realistic detector solid angle coverages and energy thresholds, was matched to 
the observed intensity of the coincidences in opposite detectors, and the intensity 
of each channel was determined by the total initial intensity of the simulation 
before applying the experimental restrictions or the experimental conditions. 
The branching ratios were obtained using for normalization the /3a activity 
from the decay of the daughter nucleus ^^Be (Ti/2=13.81(8)s (1^), measured 
without any time restrictions. The branching ratio of the ^^Be (3a channel 
following ^^Li beta decay is obtained from the product of two terms, the ^^Li (3- 
feeding to the first excited state in ^^Bc assuming that the /3- feeding to the ^^Be 
ground state is negligible, and the ^^Be /3-delayed branching ratio of the ^^Be 
(3a channel. The first term, the feeding to the ^^Be 320 keV state is taken as 



the weighted average of the values from [9|, |28|, |29|, |30|, |3l| , resulting in 7.4(3)%. 
The second term, the ^^Be (3a channel has only been determined once before 
(l9l |. with a value of 2.9(4)%. The resulting branching ratios are shown in the 
last column of Table [TJ The sum of the branching ratios to the 16 and 18 MeV 
states, 0.31(5)%, is comparable to the previously pubhshed values of /3-feeding 
to the 18 MeV state in "Be, 0.30(5)% [H and 0.39(7)% ^. 

The B(GT) values are calculated from the branching ratios. This approach 
can be problematic for broad levels, as the B(GT) is normally defined as a factor 



in the R-matrix formalism (see e.g. |2J, [33])- Instead, we used a modified ft- 



value, where the Fermi function was substituted by its average over the state 
resonance shape, which was modeled using the single-level single-channel R- 



matrix formalism 2J]. In the case of the 18 MeV state in ^^Be we included 
the neutron emission channels to states in ^''Be, with branching ratios taken 
from [HI, and the /3t value obtained in this experiment [lOj. The resulting 
B(GT) values are shown in Fig. [3I^b), represented as solid bins. The B(GT) 
values for levels below the charged particle emission thresholds are included for 
completion, and represented as open bins. These B(GT) values were calculated 
from the branching ratios of ll|, llO| , assuming the levels are narrow and using 
the regular /i-value. The filled part of the 10 MeV bin in Fig. [3^b) shows 
the B(GT) deduced from the branching ratio corresponding to charged particle 
emission from the 10.6 MeV state, taken from 

Fig. [Sja) shows the calculated B(GT) strength distribution to excited states 
in ^^Be from shell model done using unrestricted excitation space in the p and 
sd shells (for details see [1^ H^)- The inset panels in Fig. [3l^a) shows the 
partial B(GT) distribution assuming feeding only to 3/2~ or 5/2~ states in 
^^Be, respectively. The shape of the calculated B(GT) distribution to 3/2~ 
states at high excitation energy in ^^Bc (E>10.6 MeV) follows the trend of the 
experimental B(GT) values, supporting the 3/2~ spin and parity assignment 
for these states given in this work. The experimental B(GT) strength for the 
16 MeV state in ^^Be is low compared to our shell model prediction and to 
the B(GT) strength calculated in [li|. It should be stressed that the B(GT) 
values for the ^^Be excited levels given in this work correspond to decays were 
at least two charged particles are involved in the breakup. However, one should 
be aware that there is previous evidence of a feeding of up to 6% to levels in the 
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region between 14-17 MeV excitation energy in ^^Be decaying through neutron 
emission to low excited states in ^"Be [ll|. This could explain the difference 
between the expected B(GT) distribution in this region and the experimental 
one. Unfortunately, with our setup, we just enhance the two-body and three- 
body decay channels involving at least two charged particles. 

The study of the ^^Be*^a-|-^He^n-|-a-|-^He breakup channel presented in 
this Letter shows the power of the full kinematic analysis when three body 
channels participate in the decay. The presence of five body decay channels can 
obscure the signature of the three body ones, unless one has a suitable setup 
as the one used in this work. Full kinematics was achieved by detecting the a 
and the ^He in coincidence, and reconstructing the neutron energy and emission 
angle from energy and momentum conservation. This allowed us to disentangle 
the presence of this three body channel in the decay of two excited states in 
^-^Be, a known state at 18.4(3) MeV, and a state at 16.3(1) MeV, previously 
unreported in the literature. Although no experimental evidence of such state 
in ^^Be existed, previous shell model calculations 1^, 33| predicted /3~ strength 
in this region. The study of the angular correlations in the breakup of these two 
states indicates a 3/2~ spin assignment for both states. 
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Table 1: Level centroid and reduced widths used in the R-Matrix description of the states 
modeled in the Monte-Carlo code. The width F was obtained from a gaussian fit of the 
R-matrix peak directly. The last column lists the branching ratios determined in this work 
following ^^Li /3-decay. The ^^Li activity was deduced from the branching of the /3q decay 
channel of the daughter ^^Be. 
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Figure 1: (color online) Left: (a) scatter plot of the individual a, ^He and (reconstructed) 
neutron energies plotted versus the ^^Be excitation energy, recorded within the first 60 ms 
after the proton impact. The charged particles were measured in opposite detectors, and the 
energy of the neutron reconstructed assuming 3-body breakup. The line of 10 MeV offset and 
11/10 slope corresponds to the neutron emission of the ^^Be 10.6 MeV state @,[l3. The two 
dotted lines indicate the gate around the 11/7 interval corresponding to breakup through the 
^He(gs) resonance. Right: (b) reconstructed ^^Be excitation energy spectrum corresponding 
to the events inside the interval shown in the left panel. The good fit obtained in the Monte- 
Carlo simulation of the breakup process is shown in red dashed line over the data. The 
''^He+a break up contribution of the ^^Be states at 16 and 18 MeV is shown with hatched 
area histograms, and the black dash-dotted (bellow the hatched area) and green dotted lines 
correspond to the 2a3n breakup of the for ^^Be 16 and 18 MeV states respectively . 



11 








50 100 150 

e (deg) 








50 100 150 

e (deg) 



Figure 2: (Color online) Left: (a) the angular distribution for events gated on the "Be 
18 MeV state is shown. The Monte-Carlo simulation of the most favored spin and parity 
assignments is overlaid. The red dashed histogram corresponds to spin J=3/2 for the ^^Be 18 
MeV state, which is the favored spin value from the test. Right: (b) angular distribution 
for events gated on the ^^Be 16 MeV state. Again, the most favored Monte-Carlo simulation 
is shown, in this case J=3/2 (red dashed histogram) for both states at 16 and 18 MeV in 
"Be. 
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Figure 3: Left: (a) shows the B(GT) distribution to excited states in ^'^Be from a shell model 
calculation using the unrestricted excitation space in the p and sd shells [33l | . The upper panel 
shows the B(GT) distribution for 3/2~ states in ^^Be. The lower panel shows the distribution 
for 5/2~ states in ^'^Be. Right: (b) shows the B(GT) transition values to excited states in 
^^Be is shown. The values obtained in this work are shown in solid bins. The values to 
states below the charged particle emission thresholds are calculated from the branching ratios 
reported in [lll[lol| and assuming the levels are narrow and shown in open bins. The overall 
shape of the calculation for 3/2~ states agrees well with the experimental distribution, further 
supporting the 3/2 ~ assignment discussed from the angular correlations. 
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